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VARADHAN ESTIMATES FOR ROUGH DIFFERENTIAL EQUATIONS 
DRIVEN BY FRACTIONAL BROWNIAN MOTIONS 
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Abstract. In this work we study rough differential equations driven by a fractional 
Brownian motion with Hurst parameter H > j and establish Varadhan's small time 



Oj' estimates for the density of solutions of such equations, under Hormander's type condi- 



tions 
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cn ■ 1. Introduction 

Let B = {B^,...,B'^) be a d-dimensional fractional Brownian motion with Hurst pa- 
, , , rameter H > ^, that is, i? is a M"^- valued centered Gaussian process with covariance 

_c^.: HBlBi) = ^(t^^ + .2^ - It - s|2^)5,,-. 

A straightforward application of the Kolmogorov continuity theorem shows that the Hurst 
parameter H controls the path regularity of B: the sample paths of B are almost surely 
locally 7-Holder continuous for all 7 < if. When H = ^, B is a. standard Brownian 
motion. 

We are interested here in the following family of stochastic differential equations driven 
by 5: 

XI = X + eS^ I V^iXl)dBl, £G(0,1), (1) 

where Vi, . . . , V^ are C°°-bounded vector fields on M". When H > ^, the above equation is 
understood in the sense of Young's integration. In that case, existence and uniqueness of 
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solutions are well-established, for instance, in |23] and j2^. When ^ < if < |, equation 
([1]) is interpreted in the framework of rough path theory (see [T0| [20]). Existence and 
uniqueness of solutions in this case can be found, for example, in [20]. In particular, 
when H = ^, this notion of solution coincides with the solution of the corresponding 
Stratonovitch stochastic differential equation. 

Once equation ([1]) is properly interpreted and solved with a unique solution, a natural 
question to address and most relevant to our purpose is the existence of a (smooth) density 
function of the random variable X^. In the regular case H > ^, the authors proved (see [T]) 
that under the classical Hormander's condition on the vector fields, the law of X^ admits 
a smooth density with respect to the Lebesgue measure. The existence and smoothness 
of density function is more involved in the rough case. Still under Hormander condition, 
the existence of a density function for ^ < iJ < ^ is due to [7]. The smoothness of this 
density is proved in [12] for if > |, conditioned on the integrability of the Jacobian of the 
system which is established later in [9]. Finally, smoothness of the density for all if > | 
is proved in [8]. 

With the existence of the density in hands, the next step is to study some of its basic 
properties. Small time asymptotics in the regular case H > 1/2 have been studied in 
[2]. Gaussian or sub-Gaussian upper bounds have been established in [3] and j^. The 
positivity of the density is proved in [1]. In this work, we are interested in Varadhan type 
estimates for this density. 

Throughout our discussion, we will assume some uniform hypoellipticity condition in- 
spired by Kusuoka [T7] (see Hypothesis 13.41 for details). The main problem we are in- 
terested in is to establish a version of Varadhan's estimates for the density of Xf. More 
precisely, introduce the following functions on R", 

^^^y^ = ^^r^^ oll^ll.^^' ^^^ dl{y)= inf -\\h\\%, 

'S>i {h)=y I $1 (/i)=y,det 7,5^ (,,, >0 Z 

where $ is obtained by solving the ordinary diferential equation driven by Cameron- 
Martin paths: 






i=l>^0 



Our main result is the following. 
Theorem 1.1. Let us denote by Pe{y) the density o/Xf. Then 



and 



liminf£MogPe(y) > -d%{y), 



\iui sup e'^ log Pe{y) < —d^{y). 



Moreover, z/inf<i,^(/,)=y,det7<i,^(j,)>o det 7<j,,(/,) > 0, th 



en 



\im.e'^\ogPe{y) = -d\{y). 

e\.Q 
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In the classical case when i? is a standard Brownian motion, these results were studied 
by numerous authors, including Leandre [18], Ben Arous and Leandre [6] and Leandre and 
Russo [in]. Our result is obviously an extension of the classical result to the fractional 
Brownian motion case, in which many recent developments in rough paths theory are 
employed. In particular, we would like to mention the recent breakthrough [9] in the 
integrality of the Jacobian of equation ([1]). This result leads to the integrability of the 
Malliavin derivative JDXf and inverse Malliavin matrix 7^^ of X^. One of the main 
technical difficulties of this work is to show that under the uniform hypoelliptic condition, 
for any fixed r > 

for some constant Cr depending on r and constant / depending on the structure of vector 
fields V^j's which is given in Hypothesis 13.41 

Finally, let us mention that some small-time asysmptotics of density function results 
have already been studied in the smooth case H > ^ ioi elliptic systems in [2] and |15| . 
These results clearly imply the Varadhan's asymptotics stated in Theorem 11.11 In the 
rough case | < if < |, the Laplace approximation is obtain by Inahama |14| . which is 
along the same line of research as the current paper. 

The paper is organized as follows. In section 2, we give the necessary background on 
rough path and Malliavin calculus that will be needed throughout the paper. We state 
and prove our main result in Section 3. 

2. Preliminary material 

For some fixed if > |, we consider (fi, J-", P) the canonical probability space associated 
with the fractional Brownian motion (in short fBm) with Hurst parameter H. That is, 
Q = Co([0, 1]) is the Banach space of continuous functions vanishing at equipped with 
the supremum norm, J-" is the Borel sigma-algebra and P is the unique probability measure 
on Q such that the canonical process B = {Bt = [B], . . . , Bf), t G [0, 1]} is a fractional 
Brownian motion with Hurst parameter H . In this context, let us recall that i? is a 
d-dimensional centered Gaussian process, whose covariance structure is induced by 

R {t, s) := E [Bi B{] = - [s^^ + f^ - \t - s^^) , s, t G [0, 1] and j = 1, . . . , d. (2) 

In particular it can be shown, by a standard application of Kolmogorov's criterion, that 
B admits a continuous version whose paths are 7-Holder continuous for any 7 < if . 

2.1. Rough path. For A^ G N, recall that the truncated algebra T^{W^) is defined by 

N 
r?i=0 

with the convention (R'^)®° = M. The set T^{W^) is equipped with a straightforward 
vector space structure, plus an operation (g> defined by 

N 

Timig ^h) = J2 -^m-kig) ® TTk{h), g,he T^(M'^), 

fc=0 
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where vr^ designates the projection on the mth tensor leveL Then (T^(M"'), +, ®) is an 
associative algebra with unit element 1 G (R'^)'^°. 

For s < t and m > 2, consider the simplex A^ = {(mi, . . . ,'Um) ^ [-s,^]*"; Mi < ■ ■ • < 
Um}, while the simplices over [0, 1] will be denoted by A"*. A continuous map x : A^ — >■ 
T^{W^) is called a multiplicative functional if for s < m < t one has x^^^ = x^^u ® x^^f. An 
important example arises from considering paths x with finite variation: for < s < t we 
set 

<t= Yl (f^dx''---dx'Aei,^---®ei^, (3) 

l<Ji,...,i™<d ^"^^" ^ 

where {ei, . . . , e^} denotes the canonical basis of W^, and then define the truncated sig- 
nature of X as 

N 

Sn{x) : A^ ^ r^(R'^), {s,t) ^ SN{x)s,t := 1 + J2<t- 

m=l 

The function Sn{x) for a smooth function x will be our typical example of multiplicative 
functional. Let us stress the fact that those elements take values in the strict subset 
QN^^d^ C T^iM."^) given by the group-like elements 

G'^(M^)=exp®(L^(R'^)), 

where L^{M.'^) is the linear span of all elements that can be written as a commutator of 
the type a^b — bi^ a for two elements in T^iM"^). This set is called free nilpotent group 
of step A^, and is equipped with the classical Carnot-Caratheodory norm which we simply 
denote by | ■ |. For a path x G C([0, 1], (j'^(R"')), the p-variation norm of x is defined to be 

^ i/p 

|x||p_var;[0,l] = SUp | V Ix^T^ ® X^^^J^ 

nc[o,i] W 

where the supremum is taken over all subdivisions 11 of [0, 1]. 

With these notions in hand, let us briefly define what we mean by geometric rough path 
(we refer to [TU| I2U] for a complete overview): for p > 1, an element x : [0, 1] — )■ G^pJ (R*^) 
is said to be a geometric rough path if it is the p-var limit of a sequence S^p^ (x"*) of lifts 
of smooth functions x"^. In particular, it is an element of the space 

^p-var;[0,l]Q0,l],GLPj(R'^)) = {x G C([0, 1] , G^^J (R'^)) : ||x||p„.ar;[0,l] < Oo}- 

Let X be a geometric p- rough path with its approximating sequence x'", that is, x™ is a 
sequence of smooth functions such that x™ = S^pj (x"^) converges to x in the p-var norm. 
Fix any 1 < g < p so that p~^ + q~^ > 1 and pick any h G C'^~™''([0, 1],R"'). One can 
define the translation of x by h, denoted by T/i(x) by 

T^(x)= limSipiix'^ + h). 

n— >-oo 

It can be shown that T/i(x) is an element in C^~™''([0, 1],G^'''^(M.'^)). Moreover, one has 
the following continuity property. 

Lemma 2.1. For any 1 < q < p so thatp~^+q^^ > 1, let{^,h) G Cp~^^''([0, 1], GLpJ(K'^))x 
(-g-varQQ^ 1],R'^). Denoted by Tft(x) G Cp-™''([0, l],G\-P\W^))the translation ofx by h. We 
have 
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(1) There is some constant C depending only on p and q, 

\\Th{y^)\\p-va.T;[0,l] < C'(||x||p_var;[0,l] + || ^||g-var;[0,l])- 

(2) The rough path translation (x, h) i— !■ T/j(x) as a map from 

(^p-varQQ^^j^^LpJ(Rrf)) X C^-™>^([0,1],M'^) ^CP~"^'^([0,1],GLpJ(R'^)) 

is uniformly continuous on bounded sets. 

Remark 2.2. A typical situation of the above translation of x by /i in the present paper is 
when X = B, the fractional Brownian motion lifted as a rough path, and /i is a Cameron- 
Martin element of B. In this case, we simply denote Th(B) = B + h. 

According to the considerations above, in order to prove that a lift of a d-dimensional 
fBm as a geometric rough path exists it is sufficient to build enough iterated integrals 
of i? by a limiting procedure. Towards this aim, a lot of the information concerning B 
is encoded in the rectangular increments of the covariance function R (defined by ([2])), 
which are given by 

Rt = E [{Bl - Bl) {Bl - Bl) 
We then call 2-dimensional p-variation of R the quantity 

i/p 



V, 



(RY ^snp I \T\RnX\\j 



[s,),{tj)En 



where 11 stands again for the set of partitions of [0, 1]. The following result is now well 
known for fractional Brownian motion: 

Proposition 2.3. For a fractional Brownian motion with Hurst parameter H, we have 
Vp{R) < oo for all p > 1/{2H). Consequently, for H > 1/4 the process B admits a lift B 
as a geometric rough path of order p for any p > 1/H . 

Proof. The fact that Vp{R) < oo for all p > 1/{2H) is the content of [TU| Proposition 
15.5]. The implication on the rough path construction can also be found in |im Chapter 
15]. 

n 

2.2. Malliavin Calculus. We introduce the basic framework of Malliavin calculus in this 
subsection. The reader is invited to read the corresponding chapters in [21] for further 
details. Let S be the space of M"^- valued step functions on [0, 1], and Ti the closure of S 
for the scalar product: 

d 

((l[0,ti]'"' ■ 'l[0,id])'(l[0,si]5--- ,l[0,Sd]))-H = }^R{ti,Si). 

We denote by K^ the isometry between T-L and L'^{[0, 1]). When H > ^ it can be shown 
that L^/-f^([0, 1],M^) C H, and when I < H < ^ one has 

C^ cncL\[0,l]) 

for all 7 > I - i7. 
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We remark that l-i is the reproducing kernel Hilbert space for B. Let J^ be the 
Cameron-Martin space of -B, one proves that the operator TZ := TZh '■ 'H — ?■ J^ given by 

7^^:= [ KH{;s)[K*Hij]{s)ds (4) 

defines an isometry between "H and J^. Let us now quote from [T0| Chapter 15] a result 
relating the 2-d regularity of R and the regularity of J^. 

Proposition 2.4. Let B be a fBm with Hurst parameter ^ < H < ^. Then one has 
Jif C C"^'^^ for p > {H + 1/2)"^. Furthermore, the following quantitative bound holds: 



{v,{R)y/^- 

Remark 2.5. The above proposition shows that for fBm we have J^ C C"™'' for p > 
{H + 1/2) ~^. Hence an integral of the form J h dB can be interpreted in the Young sense 
by means of p- variation techniques. 

A J-'-measurable real valued random variable F is then said to be cylindrical if it can 

be written, for a given n > 1, as 

F = f (5(01), . . . ,5(0")) = /(I {<Pl,dBs), ...J (0r,rfi?.)) , 

where (f)"^ ^Ti and / : R" — t- M is a C°° bounded function with bounded derivatives. The 
set of cylindrical random variables is denoted S. 

The Malliavin derivative is defined as follows: for F E S, the derivative of F is the M*^ 
valued stochastic process (DtF)o<t<i given by 

D,F = f^0^(t)|^(5(0i),...,5(0")). 

More generally, we can introduce iterated derivatives. li F E S, we set 

ti,...,tfe ti i^fe 

For any p > 1, it can be checked that the operator D'^ is closable from S into LP(r2; T-C^^). 
We denote by D''"'^ the closure of the class of cylindrical random variables with respect to 

the norm 

1 

and 

D°° = Pi Pi d^^'P. 
p>ifc>i 

Definition 2.6. Let F = {F^, . . . , F"-) be a random vector whose components are in D°°. 
Define the Malliavin matrix of F by 

7F = ((DF\DF^)^)i<„-<„. 

Then F is called non- degenerate if'^p is invertible a.s. and 
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It is a classical result that the law of a non-degenerate random vector F = {F^, . . . , F") 
admits a smooth density with respect to the Lebesgue measure on M". Furthermore, the 
following integration by parts formula allows to get more quantitative estimates: 

Proposition 2.7. Let F = (F""^, ..., F") be a non- degenerate random vector whose com- 
ponents are in 3°°, and '-/p the Malliavin matrix of F. Let G G ©°° and (f he a function 
in the space C°°(M"). Then for any multi-index a G {1, 2, ...,?t,}'^, /c > 1, there exists an 
element Ha G ©°° such that 

ndMF)G]=nv{F)Ha]. 

Moreover, the elements H^ are recursively given by 

d 

H^,) = Y,HGiiF'r^F^) 

Ha = H^ak)iH^ai,...,ak-i)), 

and for l<p<q<oowe have 

II tr II ^ r^ ll/^,~l"p» i?ii^' wr^w 

where - = - ^-. 

p q r 

Remark 2.8. By the estimates for H^ above, one can conclude that there exist constants 
/3, 7 > 1 and integers m, r such that 

Remark 2.9. In what follows, we use Ha{F, G) to emphasize its dependence on F and G. 

2.3. Differential equations driven by fractional Brownian motions. Let fi be a 

d-dimensional fractional Brownian motion with Hurst parameter H > j. Fix a small 
parameter e G (0, 1], and consider the solution X^ to the stochastic differential equation 



d rt 



Xt = x + e^ V.{X!)dBl (5) 

where the vector fields Vi, . . . ,Vd are C°°-bounded vector fields on M". 

Proposition 12.31 ensures the existence of a lift of i? as a geometrical rough path. The 
general rough paths theory (see e.g. |10[[TT]) allows thus to state the following proposition: 

Proposition 2.10. Consider equation ^ driven by a d-dimensional fBm B with Hurst 
parameter H > ^, and assume that the vector fields V^s are G"^ -bounded. Then 

(i) For each e G (0, 1], equation ^ admits a unique finite p-var continuous solution X^ 
in the rough paths sense, for any p > jj. 

(ii) For any A > and 6 < - we have 



p 



E 



exp A ( sup \Xf\ 
te[o,i],ee(o,i] 



< oo. (6) 
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Once equation (jS]) is solved, the vector X^ is a typical example of random variable which 
can be differentiated in the Malliavin sense. We shall express this Malliavin derivative in 
terms of the Jacobian J^ of the equation, which is defined by the relation 

Setting DVj for the Jacobian of Vj seen as a function from M" to M", let us recall that J^ 
is the unique solution to the linear equation 

Jt = ld^ + eJ2 [ DV,{Xl) r, dBl (7) 

and that the following results hold true (see [7] and [23] for further details): 

Proposition 2.11. Let X^ be the solution to equation ^ and suppose the Vi's are C°°- 
bounded. Then for every i = 1, . . . ,n, t > 0, and x ^W^, we have X^'^ G D°° and 

JDiX! = 3l,VjiX!), j = l,...,rf, 0<.<t, 
where D^X^^'* is the j-th component o/D^Xj^'*, Jf = dxX^ and J^^ = Jf(J^)~^. 

Let us now quote the recent result [9], which gives a useful estimate for moments of the 
Jacobian of rough differential equations driven by Gaussian processes. 



Proposition 2.12. Consider a fractional Brownian motion B with Hurst parameter H > 
J and p > jj. Then for any rj > 1, there 
J^ defined at Proposition \2. 11\ satisfies: 



i and p > jj. Then for any rj > 1, there exists a finite constant Cr, such that the Jacobian 



E 



sup II J ||p_var;[0,l] 
ee[0,l] '^ ' ' 



c,. (8) 



In the sequel, we also need the following restatement of [T2l Proposition 4.4]. 

Proposition 2.13. Fix k eN and let {hi, . . . , hk} be any family of elements in Jif. Then 
the directional derivative Y)hi ■ ■ ■ ^hk^t exists for any t G [0, 1] . Moreover, there exists a 
collection of finite indexing sets 

{K(i,,...,i,) : (ii,...,Zfc) G {1,...,4^|, 

such that for every j G {l,..,n} we have 

Dh. . . . ^h.X^t'' (9) 

= E E / fr{ti)...fntjnTiit)dhY(ti)...dh',^{tk), 

for some functions f^'"^ which are in C'P'^^"^ ([0, 1] , M) for every p > jj, i and m, i.e. 
U(.„...,,),|i,...,,}^ U„^K(n,...,.) ^f^'"^ ■.i=l,:,k+l}C C^— ([0, 1] , M) . 

Furthermore, there exist constants C > and a > 1, which depend only on m such that 
for all t G (0, 1] 
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for every i = 1, . . . ,k + 1, every m G K(i^.,,.^i^) and every {ii, . . . ,ik) G {1, . . . ,d} , where 
we have set M" = (X% P, (P)"^). 

Remark 2.14. The fi'"^^s in the above proposition are constructed recursively from d^V{X^), 
Jl, {Jl)~^ and niultiphcations of them through iterated integrals with respect to Bg. By 
the discussion in [TOl Chapter 11], components of d^V{X^), 3^ and (J^)~^ are smooth with 
respect to e in C^"'"'"'([0, 1],M). Moreover, derivatives of dl^V{X^),Jl and (J^)""^ (with 
respect to e) satisfy linear equations similar to ([Tj) and by the same technique used in 
proving Proposition 12. 121 one can show that the derivatives are in L^{F) for all r > 1 and 
uniform in e G (0,1]. Hence /;^'™ is differentiable with respect to s in C^~™''([0, 1],]R). 
Moreover for all r > 1, one has 



df! 



E sup 

ee(o,i] 



< CX) 

p—var;[0,l] 



d£ 
for all / and m. 

Let $ : J^ — 7- C([0, 1],]R") be given by solving the ordinary diferential equation 

^,{h)=x + Yl / Vii^sih))dhi. (11) 

Following Proposition 12. 4^ Proposition 12.131 and Remark 12. 14^ we have 
Proposition 2.15. For each h G ^ , one has 

lim - ($i(efi + /i) - ^i{h)) = Z{h), (12) 

e4,0 £ 

in the topology of 3°°. 

Proof. We need to show that the convergence stated in the proposition takes place in 
II ■ ||fc,r for any A; > and r > 1. First note that Yf := ^ti^B + h) satisfies equation 

Y,^ = x + J2 [ V,{Y^)d{eBl + h^). (13) 

When H > ^, note that for any integer k > and fixed Si, ..., s^, D^^ ^ Y^ satisfies an 
linear equation. By the discussion in [T0| Chapter 11], D^ ^ Y^ is differentiable with 
respect to e as a random vector. The fact that D^^ ^ Y/ is also differentiable with respect 
to the norm E|| ■ ||^»fe follows from the fact that || ■ ||^(g>fe is controlled by the sup- norm 
and the integrability of the system. 

Next, we focus on the case H < ^, and divide the proof into three steps. 
Step 1: When k = 0. It is known that Y"^ is smooth (path-wise) with respect to e in the 
topology of CP-''"'^([0, 1],R'') (cf. ^ Chapter 11]). To see that the convergence in (^^ 
also take place in L^iF), one only needs to note that 

e e Jq d9 

where the above is considered to be an equation in C^~'""''([0, 1]; M"). Since dYf /d6 satisfies 
a linear equation, \\dY-l^ /d6\\p^,uar;[o,i] is uniformly integrable with respect to ^ G [0, 1] in 
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U'{F) for all r > 1. Hence one can conclude that \{Yf — Y^)/e\ is uniformly integrable in 
L^(F). The claimed convergence follows for the case A; = 0. 

Step 2: When k = I. Denote by Jf = -j^, the Jacobean of Yf. We have for any hi G J^ 



^k^Yi 



dx 

e l' JKJ:r'V,{Y:)dh\{s). 

'0 



Let /J'^ = sJl{Jg) ^Vi{Yg), the integrand in the above integral. It can be shown that /^ 
is smooth with respect to e in C^~'""'^{[0, 1]; M"). In particular, we have 



£,1 



de 



\e=0 



1 fe,l 



n^-f. 



0,1 



df^ 



s 



<c 



<c 



^ 

fe,l _ fO, 

Jo Jo 



dn 



de 1^=0 



£,1 



p£,l _ fO,! 

/o /o 



de ''^=0 






c?/. 



e,l 



dfo 



J S J ! 



0,1 






p— var;[0,l] 



p— ■uar;[0,l] 



|"-l||i}-mr;[0,l] 

I^iIIjT' — > as £ 4- 0. 



This implies that the convergence 

D,,F/ - D.^F^o rfD,, Fi- 



de 



le=0 



^^ as £ i 0, 



is uniform in hi G Jif. Hence 



BYf - Dyo dBY' 

-^ 



de 1^=0 



in T-L almost surely. The fact that the convergence is also in || ■ || i,r- follows from the uniform 



integrability of 



pyf-py" 



dDY-, 



de 1^=0 



m e. 



H 



Step 3: Now we proof for general k > 1. Note that for hi, . . . ,hk G Jif, the directional 
derivative Y)hi ■ ■ ■ D/ij.F/ exists for any t G [0, 1] . Moreover, by Proposition 12.131 (with a 
slight modification) there exists a collection of finite indexing sets 

|K(i,,...,i,) : {ii,...,ik) e {l,...,d}^j , 
such that for every j G {1, .., n} we have 



Dj,...Dft.ri''' 



E E 

il,...,ife=l mSK 



(14) 



£,m 



0<ti<---<t..<l 



/['^ (to . . . f^ (t,) fiz (1) dw,^{ti) . . . dhi^it,), 



for some functions /^^'™ which are in ^"^"^([0,1],]^) for every i and m. By Remark 
12.141 each fj'"^ is smooth with respect to e in C^~'"'^''([0, 1]; M) with uniform integrable 
derivatives. Now by a similar argument to that in Step 2 the proof is completed. D 
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Finally, we close the discussion of this section by the following large deviation principle 
that will be needed later. 

Theorem 2.16. Let $ be given in fTl\} . which is a differentiable mapping from M" to 
C([0, 1],M"). Denote &?/ 7$i(/i) the deterministic Malliavin matrix of^i{h), i.e., 7^ z^-, = 
(D$'^(/i), D${(/i))'^, and introduce the following functions on M" and M" x M, respectively 

^'^y^^Jjlt oll^ll.^' and J^(y, a) = inf -|l^ll^- 

Recall that Xf is the solution to equation ^ and 7xe is the Malliavin matrix of Xf. Then 
(1) XI satisfies a large deviation principle with rate function I{y). 
(2)The couple {Xf,'-fx^) satisfies a large deviation principle with rate function Iii{y,a). 

Proof. Fix any p > -^. It is known (see [10]) that B as a G'LpJ(]R'')-valued rough path 
satisfies a large deviation principle in p-variation topology with good rate function given 
by 

kWhP^ iihejr 



J(h) = i 2\\"'\\je 

1+00 otherwise. 

It is clear $i(-) : GLpJ(M'=') -^ W is continues. 

Moreover, by a similar argument as in Proposition I2.15[ one have for aA\ H > ^ 

D.$i(-) :G'LpJ(M'^)^H 

is continues. Hence 7(j>i(.) : GLpJ^M'^) — t- M is continues for all H > ^. 

Now note that Xf = $i(eB) and (X^,7xe) = ($(£:B),7$^(eB)), the claimed result 
follows from the contraction principle. D 

3. Varadhan Estimates 

Recall that we are interested in a family of stochastic differential equations driven by 
fractional Brownian motions B (with Hurst parameter H > ^) oi the following form 

Xt = x + eY,i V,iX!)dBi. 
We have defined the map $ : Jif — ?■ C[0, 1] by solving the ordinary deferential equation 

$,(/,) = x + v / v,{<^,ih))dhi. 

Clearly, we have Xf = ^^(eB). Introduce the following functions on M", which depends 
on $ 

d'{y)=Iiy)= ini hh\\%, and dl{y) = inf hh\\%. 

Assume Hypothesis 13.21 or Hypothesis 13.41 below, our main result is the following. 
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Theorem 3.1. Let us denote by Pi;{y) the density ofXf. Then 



and 



limmie^logPs{y)> -d%{y), (15) 

e4-0 



liuisup e'^ log p e (y) < —d'^{y). (16) 



Moreover, if mi^^(^h)=y,det-,i,^(y^>odet-f^^(^h) > 0, then 

lime^hgpsiy) = -d%{y). (17) 

£4-0 

A key ingredient in proving Theorem 13.11 is an estimate for the Malhavin derivative 
DXf and Malhavin matrix '-fx^ of Xf. Since it is more involved when the vector fields 
Vi's form a hypoelliptic system, we divide the estimates of theses two quantities into two 
parts: (1) when V^j's are uniformly elliptic and (2) ^j's are uniformly hypoelliptic (see 
Hypothesis | 



3.1. Elliptic Case. Throughout our discussion in this subsection, we assume that Vi, ...,Vd 
form a uniformly elliptic system. 

Hypothesis 3.2. The vector fields Vi, . . . ,Vd of equation ([5]) is uniformly elliptic, that is 

v*V{x)V*{x)v>X\v\'^, for allv.xeW, (18) 

where we have set V = {Vj)i=i^,,,^n;j=i....d o-nd where A designates a positive constant. 

Our main technical result in this subsection is the following. 

Lemma 3.3. Assume Hypothesis \3.B . For H > j, we have 

(1) supgg/Q 11 ll^ilU,r < oo for each k > 1 and r > 1. 

(2) ||7xe||r < CrS~'^ for any r > 1. 

Proof. We follow the idea in |1]. First recall that the Malhavin derivative D^Xf can be 
expressed as 

where J*^ is the Jacobian process defined by Jf = -^. We divide the proof into two cases: 
H > I and j < H < ^. 

When iJ > i, it is clear that (see [13], for example) 

sup ||X^||fcr < OO 
£6(0,1] 

for each k > 1 and r > 1. 

In the following we prove the desired bounds for the Malhavin matrix. Let 

Jo Jo 
Our bound for the Malhavin matrix 7^^ is now reduced to prove that 

y*T,y>M,\y\^ for y G M", (19) 
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for a given random variable M^ whose inverse admits moments of any order uniformly in 
e G [0, 1]. To this aim, notice first that 

y*T,y= f [ {U,f,)^,\u-v\'''~'dudv, with U^ViX'J*iiJir'ry. 
Jo Jo 

Furthermore, thanks to the interpolation inequality of jll Lemma 4.4] applied to 7 > H—^, 
we have 

/ / {fuJ.)\u-v\'''-'dudv>C^ -— ^, (20) 

^0 ^0 II/II7 

where \\f\\-y is the 7-Holder norm of / on the interval [0, 1]. As a consequence, since the 
ellipticity condition |V"(x)yp > A|yp holds true, it is readily checked that 

\M'>X\{3ir'y\'>\\\3imy\', and ||/||, < c(l + ||X^||,)(1 + || J-i||,)|y|. (21) 



Plugging these relations into ( l20|) we deduce that for every y G 

y*T;'y < c(i + ||x^||,)^(i + IK^)"!,)^!! J^l^ \y 



2 

7 IWI ) 



from which the desired result follows easily. 

Next, we prove for the case j < H < ^. We first prove claim (1) for the Malliavin 
derivatives. Thanks to Proposition 12.131 and Proposition 12.41 we have, for q > {H + |)~^ 
and hi, . . . ,hk G J^ : 

I'Dhi ■ ■ ■'Dhk^tl 

— <^l|l/fc+llloo(l + ll/r ||p-var;[0,t] • • • || /mllp-var;[0,t]) || ^1 ||g-var;[0,i] • • • || ^fc ||g-var;[0,t] 
< €2(1 + ||M^||p_var;[0,t])ll^l|lg-var;[0,t] • • • || /ifc ||g-var;[0,t] 

<c,{i + \\MTA^rm)\\hi\U'---\\hk\U- 

Hence for any r > 1, 

E||D'=X[||^«. < c (1 + E||M^||^"Xa,[o,t]) < 00. 

Next we prove the estimate for '-/x^- Let Ai';/^ = (D^X^'*, D^X^'"'). We can deduce 
that for any u G M"', 



V -fXfV 



d fl pi 

Y^ Wv^TiiXlW-u >ch I \v*T>\Xl\^ds = ch v*Mlvds. 
i^i Jo Jo 



In the above, the inequality is obtained by the fact that Ti C -^^[0, 1]. Hence 

1 r^ 

(v*jx^v)-^ < — / (v*Mlv)-^ds. (22) 

ch Jo 

Let us now derive a suitable bound for A^^: recall that the Malliavin derivative DJX[ 
can be expressed as DJX^ = £:J^i(J^)~^\^(Xj^), where J^ is the Jacobian process defined 
by Jf = ^. We have 

Mr^ = (D,xr,D,X['^-) = e2((JKJ:)-V(Xn)\ {JUJiy'ViX^. 
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Hence by the uniform integrability of J^ and (J"^) ^ in e G [0, 1], and the uniform eUipticity 
of the vector fields V^^'s, we easily bound 



E 



sup sup I — 

ee(0,l]sg[0,l] V^' 



< Cr 



for any r > 1, where A^ is the smallest eigenvalue of Ai^. Hence (122|) implies for any r > 1 
sup p| ^ ^5^' <5l<supp(— [ {v*M'sVy^ds > 6-^\ 

\v\ = l I ^ J |i)|=l [Ch Jo J 

A 

sG[0,l] \ ^ 



<¥{ sup [ch^] >6-'} <Cr,H6''. 



Now we can conclude, by |2n Lemma 2.3], that ||7xe||r < c^e ^. This yields the claimed 
result. 

D 



3.2. Hypoelliptic case. In this subsection, we extend the results in the above under a 
weaker assumption on the vector fields Vi, ■ ■ ■ , V^. We first introduce some notations. Let 
^ = {0} U Ufeli{l> 2, ■ ■ ■ , n}^ and Ai = A\ {0}. We say that J G ^ is a word of length k 
if J = (ii, ■ ■ ■ ,ik) and we write \I\ = k. If J = 0, then we denote |/| = 0. For any integer 
/ > 1, we denote by A{1) the set {I € A; \I\ < 1} and by ^i(/) the set {I G ^i; |/| < /} . 
We also define an operation * on ^ by / * J = (zi, ■ ■ ■ , ik,ji, ■ ■ ■ ,ji) for / = (zi, ■ • • ,ik) 
and J = (ji, ■ ■ ■ ,ji) in A. We define vector fields Vj/j inductively by 

V[j] = Vj ' V[uj] = [V[7] ,Vj], J = !,■■■ ,d 

Now we introduce the following uniform hypoelliptic condition, which is in force through 
out the rest of the section. 

Hypothesis 3.4. (Uniform hypoelliptic condition) The vector fields Vi,--- ,Vd are in 
C^{W^) and they form a uniform hypoelliptic system in the sense that there exist an 
integer I and a constant A > such that 

Y, {V[j]{x),u)l,.>X\\ur (23) 

holds for any x, m G M" 

Remark 3.5. It is clear that Hvpothesis 13.21 is a special case of the above Hypothesis 13.41 

The main result of this subsection is the following counterpart of Lemma 13.31 in the 
hypoelliptic case. 



Lemma 3.6. Assume Hypothesis \3.4\ For H > j, we have 



(1) sup^g(-Q ;^i ||Xf llfc^r < oo for each k > 1 and r > 1. 

(2) ||73^-e||r < CrS'"^^ for any r > 1. 
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Remark 3.7. It is clear that following the same lines in the proof of Lemma [3. 3[ one has 
the claimed estimate (1) in Theorem 13.61 Hence in what follows, we focus on establishing 
(2) of Theorem 13.61 



Under the Hypothesis 13.41 above, for any / G ^i, we can find functions Wj G C^(M' 



) ■ 



such that: 



\/[,](x)= Yl ^/(^)%](^) (24) 



holds for any x G M" 

We consider a family of SDEs indexed by e G (0, 1]: 



X^ = x + eJ2 [ V,{Xl)dBl = x + J2 [ V:{Xt)dBi 



(25) 



where the rescaled vector fields V^ are defined as V{^{x) = eVi(x). More generally, for any 
/ G Ai{l), we denote Vfj^{x) = el^lV[/](x). Note that for I G ^i(/), 

V^j^ix) =e\'\Vjix) 

= Y, el^la;/(x)y;(x) 

j£Ai{l) 

= J2 e(l'l~"l^c./(x)y[}](a;) 

j£Ai{l) 

= J2 ^f (^)^[}](^) 

JgAiH) 

where ujf^{x) = e^^^^'^-^^^ujf{x). 

It is known that for any e G (0, 1] and any t > 0, the map x — )■ X^ : R" — > M" is a flow 
of C°° diffeomorphisms (see [lO] ) . We denote the Jacobian by 

J J = dxX^. 

As being mentioned earlier, J^ and (Jf)~^ satisfies the following linear equations: 

Jt = id, + Y [ Dv;{xi)r,dBi 

i=i "^0 
and its inverse {Jl)~^ satisfy the linear equation: 

(j^')-i = id,-Y [\jT' Dv;{xi) dBi 
7=1 -^0 
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Let us introduce a linear system I3j''^{t,x) that satisfies the following linear equations: 



'ti 



(26) 



Lemma 3.8. Fix e G (0, 1]. For any I E Ai{l), we have: 

J€Ai{l) 

Proof. To sinipify the notation, let us denote 

and 

b}{t,x)= J2 f^f'{t.x)Vfj^{x). 

J€Ai{l) 

Clearly by definition, we have a]{0,x) = b'j{0,x) = Vrn(a;). Next, we show that a'j{t,x) 
and b'^{t,x) satisfy the same differential equation. Indeed, by change of variable formula, 
we have: 

da]{t,x)=i3tr\V{i^iXt)) 

d 

= Y.^-^)^^l)-'n^.v;]{xt){x)dBi 

d 

=E E -u^U^iWtr'v(j^{xi)dBi 

i=i JeAi{i) 

d 

= E E -^f:,{Xl)a^j{t.x)dBl 

i=i JeAiii) 

On the other hand, by the definition of Pj''^{t, x), we have: 
dh]{t,x)=d{ Yl Pf'\t,x)V(i,^{x)) 

KeAiil) 

= Y dpf''{t,x)Vf^^{x) 

K£Ai(l) 

=E E -^U^t) E Pj"it,x)vf^]ix)dBi 

j=i JeAiH) KeAiil) 

d 
= E E -u^UxtWj{t,x)dBl 

3 = 1 J£Ai{l) 

The result follows by the uniqueness of the differential equation. D 
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Let US introduce the following notation: for any I, J E Ai{l), we define 
M,V(t,a;) =r(l^l+l^l)^(/3^'^(-, x)l[o,](-), /3^'"(-, x)l[o,](-))«. 

In what follows, we will only consider the case t = 1 and write Mj j{x) instead of 
Mj j{l,x). We cite the paper j5l Theorem 3.5] for the following result: 

Proposition 3.9. For any p E (1, oo), 

sup E{\\{Mlj{x))j,jeMi)\\'n < oo 
ee(o,i],xeM" 

Finally, we are able to prove (2) of Lemma 13.61 

Lemma 3.10. Fix e G (0, 1] and let ^x^ be the Malliavin matrix of XI, then ^x^ is 
invertible and there exists a random variable C G L^ for p >2 such that 

C 
Remark 3.11. It follows from the above lemma that for any r > 1 

hxlWr < Cr€~^\ eG(0,l), 

for some constant Cr depending on r. 
Proof. First note that: 

= J2 (3^''it,x)JlV^j^{x) 

and 

Y^. = {D,Xl'\D,Xl^^)^. 

Hence we have the following expression for 7^=: 

7xf= Y. (/3''1-,^),/3''1-,^))«JI^[l](^)V[}](x)*(JD* 
l,JeAiil) 

Now pick u G M"', we have: 

u*^xtu= Yl (/3''1-,^),/3^'1-,^))«(^1](^),(Ji)*^)m"(^}](^),(Ji)*«)m" 
/,Je.Ai(0 

= Y M,V(x)(V^[1](x),(JD*«)m"(V^}](x),(JD*«)r" 

I,J€Ai{l) 

leAiil) 

=X^UMlAx)) Y e'\'\{V[j]ix),{Jiru)l. 
leAiil) 

>e^'KUMUx))\\\{Jiruf 
>Ae2'A„in(Mf,(x))ALn(JDhf 
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Now by choosing u the eigenvector that corresponds to AminlTxO' ^^ obtain: 

Now by the uniform integrabihty of (Jf)"^ in e G (0, 1], and proposition fl3.9p . we obtain: 

and it completes the proof. D 

3.3. Proof of Theorem 13. li Now we are in position to prove the main result of this 
paper. 

Proof of (I15j) Fix y G M". We only need to show for dj^{y) < oo, since if dj^{y) = oo the 
statement is trivial. Fix any 77 > and let h G Jif be such that $i(/i) = y, det^^(/i) > 0, 
and WhW^j^ < d'j^{y) + rj. Let / G C^{W^). By Cameron-Martin theorem for fractional 
Brownian motions, we have 

ll'»ll «f B(h) 

Consider a function x ^ C°°(]R),0 < x < l? such that x{t) = if t ^ [—2'r],2r]], and 
xit) = 1 if t G [—77, 77]. Then, if / > 0, we have 

\\h\\ M«+4r) 

E/(Xn > e S^Ex{eB{h))f{<^,{eB + h)). 



Hence, we obtain 
e^ log p 
On the other hand, we have 



sHogpM > -{^\\h\\%. + 2r])+eHogE{x{eB{h))6y{<^,{eB + h))). (27) 



E{x{eB{h))6y{<^,{eB + h))) = s^'^E ( x{eB{h))So ' ^'^'^ + ^^ *'^^^ 



Note that 

e^O e 

is a n-dimensional random vector in the first Wiener chaos with variance 7<j>i(/i) > 0. 
Hence Zi{h) is non-degenerate and by Proposition I2.15[ we obtain 



Therefore, 



lime^ \ogE{x{eB{h))6y{<^i{6B + h))) = 0. 



Letting e 4- in (127|) we obtain 

liminf£2logp,(y) > _(l||/i||2^, + 2r/) > -(4(y) + 37/). 

e4,0 Z 

Since r] > is arbitrary, this completes the proof. D 

As a direct consequence of the above proof, we have the following lower bound for the 
density function for small e G (0, 1). 
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Corollary 3.12. Assume the same conditions as Theorem \3.1[ Denote by Ps{y) the 
density of Xf. Then for all small e G (0, 1) we have 

C _4w 
Pe{y) > — e ^^, 

for some constant C > 0. 

Proof of (lisp . Fix a point y G M" and consider a function x ^ C^(]R"), < x < 1 such 
that X is equal to one in a neighborhood of y. The density of Xf at point y is given by 

My) = ^x{xt)6y{xi). 

By Proposition 12. 7[ we can write 

ExiXl)6yiXl) =E {l{xi>y}H(i,2,...,n)iXlxiXt))) 
<E|iJ(l,2,...,„)(Xl^x(^l))l 

=^{\H(l,2,...,n){Xl,x{Xl))\l{x(esuppx}) 

<F{X! G snppx)'\\Hii,..,n){Xlx{Xm\p, 
where- + - = 1. By Remark 12.81 we know that 

for some constants /3,7 > and integers k,m,r. Thus, by Lemma [3.31 we have 

hm£2log||/7(i„„,„)(X^x(^i))||p = 0. 

Finally by Theorem I2.16[ the large deviation principle for Xf ensures that for small e 
we have 

F{Xl G suppx)' < e-^^^'^f-^^— '''(^^ 
which concludes the proof. D 

Proof of (I17h . Fix a point y G M" and suppose that 

7:= inf det7$(/i)>0. 

^{h)=y,dct ■yg>{h)>0 

Let X ^ C^(]R"), < X ^ 1 be a function such that x is equal to one in a neighborhood 
of y, and g G C°°(M), < 5- < 1, such that g{u) = 1 if \u\ < ^7, and g{u) = if |m| > ^7. 
Set Gg = (7(det7xe). As before, we have 

Ex{Xl)6y{X!) = EGex{Xt)6y{X!) + E(l - Ge)xm)Sy{Xl) = h + h. 

In what follows, we estimate /i and I2 respectively. 

Estimate of Ii. Let {0n,?i > 1} be an orthonormal basis for "H. Let 

N 



i?^ = 5^i?(0„) 



i=l 



be the Karhunen-Loeve type approximation of B. Denote by B^ the lift of B^ to 
G'Lp-'(M"'). It has been shown (cf. jTU]) that B^ converges to B in the rough path topology 
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in L'^{¥) for any r > 1. Then by the continuity of the Ito's map $, we see immediately 
that Xl' = $(eB^) converges to Xf = $(£:B), as N approaches to infinity, in L^(F) 
for any r > 1. Moreover, one can show that this convergence indeed takes place in D°°, 
which can be seen by a similar (but simpler) argument to the proof of Lemma 12.151 

Now we can claim that KGsX{^i)^y{^i) — 0- Because, otherwise, based on the above 
Karhunen-Loeve type approximation and by some standard argument, one can find an 
element e/i G "H such that $(e/i) = y and < det7$(5/i) < |, and this is in contradiction 
with the definition of 7 (see [221 Proposition 4.2.1] for more details). 

Estimate 0//2: Proceding as in the proof of ( ITB]) we obtain 

E(l - GMXt)6y{Xl) =E{l{xf>y}H^i,...,n)iXl (1 - GMXm 

<E\H^^,„„^){Xl{l-G,)x{Xl))\ 



<P (xi e suppx,det7x^. > 17") ' ||if(i,...,„)(X[,x(X[; 



Finally, by Lemma 13.31 and large deviation principle stated in Theorem 12.161 for the 
couple {XI , 7xe ) , we have for any q > 1 

lim Slips'^ log pe{y) <- — inf ||/i||^ 

elO Zq <I>(fe)gsuppx,det7<s.(/i)>i7 

<-l inf 4(y). 
Zq j/Gsuppx 

The proof is completed. D 
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